Coupled quintessence with a possible transient accelerating phase 
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We discuss some cosmological consequences of a general model of coupled quintessence in which 
the phenomenological coupling between the cold dark matter and dark energy is a function of the 
cosmic scale factor e(a). This class of models presents cosmological solutions in which the Universe 
is currently dominated by an exotic component, but will eventually be dominated by cold dark 
matter in the future. This dynamical behavior is considerably different from the standard ACDM 
evolution, and may alleviate some conflicts in reconciling the idea of the dark energy-dominated 
universe with observables in String/M-theory. Finally, we investigate some observational features 
of this model and discuss some constraints on its parameters from current SNe la, BAO and CMB 
data. 



PACS numbers: 98.80.Es, 95.35.+d, 95.36.+X, 98.65.Dx 



I. INTRODUCTION 

One of the main open problems in Cosmology is to de- 
termine the physical mechanism behind the current cos- 
mic acceleration. This phenomenon has been evidenced 
by a combination of observational data [l[ and, in the 
context of the general relativity theory, can be explained 
only if we admit the existence of an exotic field, the so- 
called dark energy. The origin and nature of this exotic 
component constitute a complete mystery and represents 
one of the major challenges not only to cosmology but 
also to our understanding of fundamental physics (see, 
e.g., for more on this subject). By assuming a spa- 
tially flat geometry, this mysterious component accounts 
for (in units of the critical density) ~ 0.7 of the cosmic 
composition, a value that is of the same order of mag- 
nitude of the relative density of the cold dark matter, 
~ 0.3. However, since these dark components are usu- 
ally assumed to be independent and, therefore, scale in 
different ways, this would require an unbelievable coinci- 
dence, the so-called coincidence problem (CP). 

A phenomenological attempt at alleviating the CP is 
allowing the dark matter and dark energy to interact. 
This phenomenology in turn gave origin to the so-called 
models of coupled quintessence, which have been largely 
explored in the literature [3-6]. These scenarios are 
based on the premise that, unless some special and un- 
known symmetry in nature prevents or suppresses a non- 
minimal coupling between these components (see 0] for 
a discussion), a small interaction cannot be ruled out. 

The usual critique to coupled quintessence scenarios is 
that, in the absence of a natural guidance from funda- 
mental physics, one needs to specify a possible interact- 
ing or coupling term between the two dark components 
in order to establish a model and investigate their obser- 
vational and theoretical implications. In this concern, a 
still phenomenological but very interesting step toward 
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a more realistic interacting or coupling law was recently 
discussed in Ref. [|[ (see also [t|) in the context of models 
with vacuum decay (uj — —1). Instead of the traditional 
approach, Ref. @ deduced the new interaction law from 
a simple argument about the effect of the dark energy on 
the CDM expansion rate. Such a coupling is similar to 
the one obtained in Ref. [10( from arguments based on 
renormalization group and seems to be very general, hav- 
ing many of the previous attempts as a particular case. 

An important aspect worth emphasizing is that in the 
above analyses the interacting parameter e has been con- 
sidered constant over the cosmic evolution whereas in a 
more realistic case it must be a time-dependent quantity. 
In Ref. [ll[ the above arguments were extend for the case 
in which the interacting parameter e is a function of the 
scale factor a. The analysis of Ref. |77| . however, was 
restricted to the case in which u) = — 1, which is mathe- 
matically equivalent to dynamical A scenarios. 

In this paper, we extend the arguments of Ref. [ll| to 
a dark energy/dark matter interaction, where the dark 
energy component is described by an equation of state 
p x = ujp x with uj < 0. We explore the dynamical behav- 
ior of this class of models and find viable cosmological 
solutions for two parameterizations of e(a). In particular, 
the solutions with transient acceleration, as consequence 
of the interaction in the dark sector, are investigated in 
more detail. In order to check the observational viability 
of this general class of coupled quintessence scenarios, we 
also carry out a statistical analysis with recent observa- 
tions of type la supernovae (SNe la) along with recent 
estimates of the CMB/BAO ratio at two different red- 
shifts z = 0.20 and z = 0.35. 



II. INTERACTION IN THE DARK SECTOR 

According to current observations, the main contribu- 
tions to the total energy-momentum tensor of the cosmic 
fluid are non-relativistic matter (baryonic plus dark) and 
a negative-pressure dark energy component. By assum- 
ing that a possible interaction occurs in the dark sector, 
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FIG. 1: Evolution of the density parameters Qi (i = b,dm,x) as a function of log(a) for some selected combinations of eo 
and £, [Eq. ©] and two characteristic values of the equation-of-state parameter, w = —0.9 and w = —1.1, corresponding to 
quintessence and phantom behaviors, respectively. The case oj = —1.0 was discussed in Ref. [Tlj . 



the energy conservation equation for the two interacting 
components can be written as 
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where pdm and p Xl are the energy densities of the dark 
matter and dark energy, respectively, whereas p x is the 
dark energy pressure. As the dark components are ex- 
changing energy, dark matter density will dilute in a rate 
whose deviation from standard case, pdm oc a -3 , may be 
characterized by the function e(a), i.e., 
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where the subscript denotes current values and we have 
set ao = 1. Note that, contrarily to most analyses avail- 
able in the literature, we consider the interaction pa- 
rameter as a function of the cosmological scale factor, 
e = e(a). In what follows we also consider that the dark 
energy is described by an equation of state p x = ujp Xl 
with uj — constant < 0. 

By substituting the above evolution law into Eq. ([1}, 
we find 
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where a prime denotes derivative with respect to scale 
factor a. For u = — 1 the above equation reduces to the 



vacuum decaying scenario recently discussed in Ref. (ll| , 
whereas for lj ^ — 1 and e = const., the above expressions 
reduce to the scenario recently discussed in Refs. [J, |6j. 

We assume from now on vanishing spatial curvature 
(and neglect the radiation contribution), so that the 
Friedmann equation for this interacting dark matter-dark 
energy scenario can be written as 



n b , a a- 3 + fW a- 3+£(a) + fi*,o/(o) 
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where % = H/H a , and ri^O) ^dm,o and £l Xt o stand for, 
respectively, the current baryon, dark matter and dark 
energy density parameters. In the above equation, the 
dimensionless function f(a) takes the following form: 
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A. e(a) parameterization 



In order to proceed further and study some cosmo- 
logical consequences of the class of coupled quintessence 
scenarios discussed above, we must assume an appropri- 
ated relation for e(a). In our analysis, we consider two 
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FIG. 2: The same as in Fig. 1 for some selected values of eo, [Eq. ©]. Note that for large positive values of eo > 1.2 (Panel 
2b and 2d) , the interaction between dark matter and dark energy will drive the Universe to a new matter-dominated era in the 
future, when a — > oo. 



different parameterizations for the interacting parameter: 
e(a) =e a« =e (l + z)-«, (PI) (6) 

and 

e(a) = e exp(l-a- 1 ) = eoexp(-^), (P2) (7) 

where eo and £ may, in principle, take negative and pos- 
itive values. PI is certainly a very simple choice among 
some physically possible functional forms. Note, how- 
ever, that for negative values of £, PI blows up in the 
past, when a — > 0. Differently, P2 is a one-parameter, 
well-behaved function during the entire evolution of the 
Universe. Note also that P2 implies a weaker dark mat- 
ter/dark energy interaction in the past, as z increases. 

III. DYNAMICAL BEHAVIOR 

The time evolution of the density parameters Ob (a), 
(a) and £l x (a) can be derived by combining Eqs. ([2])- 
((5). They read: 
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where A = O dmi0 /O(, )0 , B = fi^o/fifc.o and C = 

Figure 1 shows the evolution of the density parame- 
ters as function of log(a) for PI [Eq. ((BJ]. For simplicity, 
we consider two characteristic values of the equation-of- 
state parameter, w — —0.9 and w = —1.1, corresponding 
to quintessence and phantom behaviors, respectively. In 



agreement with current WMAP results [l8j], we assume 
tt bfi = 0.0416 and £l dm ,o = 0.24. Note that, although 
currently accelerated (and, therefore, possibly in agree- 
ment with SNe la data), models with eo > and negative 
values of £ (Figs, la and Id) fail to reproduce the past 
dark matter-dominated epoch, whose existence is funda- 
mental for the structure formation process to take place. 
In both cases, the dark energy and dark matter densities 
vanish at high-z and the Universe is fully dominated by 
the baryons (for a CMB analysis in a baryon-dominatcd 
universe, see |19||). 

Regardless of the sign of eo and the values of w, well- 
behaved scenarios are obtained when £ takes positive val- 
ues (Figs, lb and le). In these cases a mix of baryons 
(< 20%) and dark matter (> 80%) dominates the past 
evolution of the Universe whereas the dark energy is al- 
ways the dominant component from a value of a < 1 
on. An interesting and completely different future cos- 
mic evolution is obtained when eo > and the parame- 
ter £ takes large positive values (> 0.8). This is shown 
in (Figs, lc and If) for £ = 1.0 and e = 0.1. Note 
that, besides having a well-behaved past evolution and 
being currently accelerating, the cosmic acceleration will 
eventually stop at some value of a » 1 (when the dark 
energy becomes sub-dominant) and the Universe will ex- 
perience a new matter-dominated era in the future, when 
a — > oo. This kind of dynamic behavior is not found in 
most coupled quintessence models discussed in the liter- 
ature, being essentially a feature of the so-called thaw- 
ing [14| and hybrid [151 ] potentials, which in turn seems to 
be in good agreement with some requirements of String 
or M theories, as discussed in Ref. [l6| (see also 



1 The argument presented in Ref. flU is that an eternally accel- 
erating universe, a rather generic feature of many quintessence 
scenarios (including the standard ACDM model), seems not to 
be in agreement with String/M-theory predictions, since it is 
endowed with a cosmological event horizon which prevents the 
construction of a conventional S-matrix describing particle inter- 
actions. 
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FIG. 3: Deceleration parameter as a function of log(a) for some selected values of eo and £ for PI and P2, respectively. Note that 
for eo > and large positive values of £ (PI) and also when eo > 1.2 (P2) the Universe will experience a new matter-dominated 
era in the future, when a — > oo. 



In Fig. 2 it is shown the same as Fig. 1 for P2 [Eq. 
0]. Note that, independently of the values of u), all 
values of eo > give rise to well-behaved scenarios in 
which the Universe had a last evolution dominated by 
a mix of baryons and dark matter and it is currently 
accelerating (dominated by dark energy) . Note also that 
for large positive values eo, e.g., e > 1.2 (Panels 2b and 
2d) the Universe will evolve to an eternal deceleration 
phase instead of the usual de Sitter phase. 

The deceleration parameter, defined as q = — aa/a 2 , is 
given by 

= 3 n b , a~ 3 + n dmfi a- 3+ ^ + (1 + Lo)n x! of{a) 
q 2 n b .oa- 3 + Oz m .oa- 3+e(a) + fi*,o/(a) 

(9) 

and shown in Fig. 3 as a function of log(a) for PI and P2, 
respectively. In agreement with our previous discussion, 
we clearly see a transient acceleration phenomenon for 
some selected values of eo and w. 

IV. OBSERVATIONAL ANALYSIS 

As we have seen, the model here discussed comprises 
a multitude of cosmological solutions for different com- 
binations of its parameters. In this section, we will 
discuss more quantitatively the observational aspects of 
this class interacting scenarios. To this end we perform 



a joint analysis involving current SNe la, CMB/BAO 
data. In our analysis, we fix fi^o = 0.0416 from WMAP 
results [18J (which is also in good agreement with the 
bounds on the baryonic component derived from primor- 
dial nucleosynthesis 20]) and consider the recent deter- 
mination of the Hubble parameter Hq = 74.2 ±4.8 [2l| in 
conjunction with the CMB constraint £ldm,oh 2 — 0.109 ± 
0.006 [18J. We use a recent SNe la compilation, the so- 
called Union sample compiled in Ref. [22[ which includes 
recent large samples from SNLS || and ESSENCE Q 
surveys, older data sets and the recently extended data 
set of distant supernovae observed with the Hubble Space 
Telescope. The total compilation amounts to 414 SNe la 
events, which was reduced to 307 data points after selec- 
tion cuts. 

Following Ref. j25j we use constraints derived from the 
product of the CMB acoustic scale 

£ A = nd A (z*)/r s (z*) , (10) 

and the measurement of the ratio of the sound horizon 
scale at the drag epoch to the BAO dilation scale, 

r s {z d )/D v (z BAO ) . (11) 

In the previous expressions, d A (z*) is the comoving 
angular-diameter distance to recombination z* = 1089 
and r s (z*) is the comoving sound horizon at photon 





FIG. 4: The results of our statistical analyses. Contours of x 2 in the plane uj — eo (left panel) and eo — f> with uj = — 1, (right 
panel) for PI. These contours are drawn for A^ 2 = 2.30 and 6.17. The best fit values are eo = —0.11 and uj = —1.04 (left 
panel), while in the case A(t)CDM (right panel) we have found eo = 0.49 and £ = 0.41. 





FIG. 5: Contours of x 2 hi the plane uj — eo (left panel) and O.dm,o — eo, with uj = — 1, (right panel) for P2. The best fit values are 
eo = —1.19 and uj = —1.16 (left panel), while in the case A(t)CDM (right panel) we have found fi<j m ,o = 0.245 and eo = 0.035. 



decoupling given by r s (z,) = f_°° j^jdz, which de- 
pends upon the speed of sound before recombination (c s ). 
Zd — 1020 is the redshift of the drag epoch (at which the 
acoustic oscillations are frozen in) and the so-called dila- 
tion scale, Dy, is given by Dy(z) = [czr 2 (z)/H(z)] 1 ^ 3 . 

By combining the ratio r s (zd = 1020)/r s (z» = 
1090) = 1.044 ± 0.019 [26] with the measurements of 
r s (zd)/Dv{z-BAo) at zbao = 0.20 and 0.35 from Ref. [27| . 
Sollerman et al. (2009) found 

fo.m = d A (z,)/D v (0.2) = 17.55 ± 0.65 

/0.35 = d A (z*)/D v (0.35) = 10.10 ± 0.38 . 
In our analysis, we minimize the function \t — XsNe + 

XCMB/BAO' Where XCJMB/BAO = ^-^ Z \ S ) ~ ^.2? /og. a + 

[fo.35(z\ s ) — /0.35] 2 /°"o.35 an( i s stands for the model pa- 
rameters. This total \t function, therefore, takes into 
account all the observational data discussed above. 

The results of our statistical analyses are shown in 
Figs. 4 and 5. We show 1 and 2a confidence regions 
in the parametric spaces: ui — eo and eo — £ f° r PI an d 
uj — eo and fldm. — ^o for P2 that arise from the joint 



analysis described above. Note that in all panels both 
negative and positive values for the interacting parameter 
are allowed by these analyses. Physically, this amounts 
to saying that not only an energy flow from dark energy 
to dark matter (eo > 0) is observationally allowed but 
also a flow from dark matter to dark energy (eo < 0) [see 
Eq. @]. In right panel of Fig. 4 we show the analysis 
for eo and £ by fixing the dark energy EoS at uj = —1, 
which is fully equivalent to the vacuum decay scenario 
proposed in Ref. 
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As expected, we note that the 
current observational bounds on £ are quite weak since 
it appears as a power of the scale factor in the energy 
density [Eqs. @ and QJJ]. When £ takes more negative 
values eo — > 0, i.e., this scenario behaves very similarly 
to the standard ACDM model. Note also that the £ldm,o 
parameter (Fig. 5 - right panel) is very well bounded 
by observations. In the case of PI, the best-fit found 
are eo = —0.11 and uj = —1.04 (left panel), whereas for 
A(t)CDM model (right panel) we have found eo = 0.49 
and £ = 0.41. For P2 we have found the following best- 
fit values eo = —1.19 and uj = —1.16 (left panel), and 
ft rfm ,o = 0.245 and e = 0.035 (right panel - A(i)CDM 
scenario). 
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V. FINAL REMARKS 

In this paper, we have investigated a general class of 
models with interaction in the dark sector whose evolu- 
tion law of the dark energy is deduced from the effect 
of the same on the CDM expansion rate. Contrary to 
most similar analyses available in literature, we consider 
a more general case in which (i) the interaction term e is 
a function of the scale factor and (ii) the EoS parameter 
may take any value (w < 0). We have shown that many 
previous phenomenological approaches discussed in the 
literature are particular cases of our approach. 

We have investigated the dynamical behaviour of this 
scenario and found a number viable of cosmological so- 
lutions for two parametrizations of e(a) (Figs. 1, 2 and 
3). In the first case (PI), when eo > and the £ param- 
eter takes large positive values (> 0.8) we have found 
solutions with transient acceleration, in which the dark 
matter-dark energy interaction will drive the Universe to 
a new matter-dominated era in the future. P2, although 
depends only on the eo parameter, also supplies solutions 



with transient acceleration, when eo > 1.2. As mentioned 
earlier, this kind of solution seems to be in agreement 
with theoretical constraints from String/M theories on 
the quintessence potential V(4>) or, equivalently, on the 
dark energy equation-of-state w, as discussed in Ref. [lf|. 

By combining recent data of SNe la (Union sample) 
with the so-called CMB/BAO ratio at two redshifts, z = 
0.2 and z = 0.35, we have shown that that both an energy 
flow from dark energy to dark matter as well as a flow 
from dark matter to dark energy are possible. For the 
two different parameterizations discussed here, we have 
also investigated A(i)CDM scenarios for which oj = — 1. 
We have also found that positive values of £ are largely 
favored over negative ones. 
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